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I. INTRODUCTION 



Treating the renormalized stress-energy tensor as the source term of the semiclassical 
Einstein field equations, one could, in principle, determine the back reaction of the quantized 
fields upon the spacetime geometry of a black hole unless the (unknown) quantum gravity 
effects become dominant. Mathematical difficulties encountered in the attempts to construct 
characteristics of the vacuum polarization in a concrete spacetime are well known and since 
the back reaction equations require knowledge of the functional dependence of stres-energy 
tensor of the quantized field, {T^) ren , on a wide class of metrics, purely analytical treatment 
of the problem is impossible. It is natural therefore that much effort have been concentrated 
on developing approximate methods. 

The vacuum polarization effects of the massive fields in the curved background has been 
studied by a number of authors [1-18]. It has been shown that for sufficiently massive fields 
(i.e. when the Compton length is much smaller than the characteristic radius of curvature, 
where the latter means any characteristic length scale of the geometry) the asymptotic 
expansion of the effective action in powers of mr 2 may be used. It is because the nonlocal 
contribution to the total effective action can be neglected, and, consequently, the vacuum 
polarization part is local and determined by the geometry of the spacetime in question. 

In the black hole spacetimes the vacuum polarization of the massive scalar, spinor, and 
vector fields have been constructed in a series of papers by Frolov and Zel'nikov in the 
vacuum type-D geometries [1-4]. They used general framework of the Schwinger - DeWitt 
method [12-18] and constructed the first order of the effective action, omitting the terms 
that do not contribute to a Ricci-flat spaces. Using a different method, Anderson, Hiscock, 
and Samuel evaluated the approximate (T^) ren of the massive scalar field with arbitrary 
curvature coupling for a general static, spherically symmetric spacetime and applied obtained 
formulas to the Reissner- Nordstrom spacetime [8]. Their approximation is equivalent to the 
Schwinger - DeWitt expansion; to obtain the lowest (i. e. mr 2 ) terms, one has to use 
sixth-order WKB approximation of the mode functions. Numerical calculations reported in 
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Ref.[8] indicate that the Schwinger-DeWitt method always provide a good approximation 
of the renormalized stress energy tensor of the massive scalar field with arbitrary curvature 
coupling as long as the mass of the field remains sufficiently large. 

The aim of this paper is to construct the renormalized stress-energy tensor of the massive 
scalar with arbitrary curvature coupling, spinor, and vector fields in the geometry of the 
Reissner-Nordstrom black hole. To our knowledge the spinor and vector fields have not 
been discussed earlier. We shall achieve this using the standard result of the theory of 
quantized massive fields in the curved background that connects the coincidence limit of the 
HDSM (Hadamard-Minakshisundaram-DeWitt-Seely) coefficient [a 3 ] with the lowest order 
of the one-loop effective action and consequently with the regularized stress-energy tensor 
[1-4, 19-23]. Indeed, functionally differentiating the effective action we obtain a general 
and rather complex expression for the renormalized stress-energy tensor that is valid in any 
spacetime. Then we specialize thus obtained formulas to the spacetimes with vanishing 
curvature scalar and apply the result to the Reissner- Nordstrom geometry. We show that 
for the scalar field the resulting (T^) ren is identical with the tensor obtained earlier by 
Anderson, Hiscock, and Samuel and that in the limit of vanishing electric charge it reduces 
to the stress-energy tensor constructed by Frolov and Zel'nikov. 

There is another important limit of the general Reissner-Nordstrom geometry that yields 
the extremal black hole. Expanding the near-horizon region of such a geometry into a 
whole manifold one obtains the Bertotti-Robinson spacetime actively studied recently. We 
construct the stress-energy tensors in the Bertotti-Robinson spacetime taking appropriate 
limits in our general formulas and analyse the conditions under which this geometry is a 
self-consistent solution of the semiclassical Einstein field equations. Analyses carried out in 
the Bertotti-Robinson spacetime yield similar results. 

The effective action approach that we employ in this paper requires the metric of the 
spacetime to be positively defined. Hence, to obtain the physical stress-energy tensors one 
has to analytically continue at the final stage of calculations their Euclidean counterparts. It 
should be stressed once again that the method, when applied to the rapidly varying or strong 

3 



gravitational fields, breaks down and that its massless limit is contaminated by nonphysical 
divergences. 



II. THE EFFECTIVE ACTION 

We begin with a short description of the method. More detailed presentation may be 
found in [3,4, 21-23]. Our notation corresponds to those of Refs. [21-23]. Consider the 
elliptic second-order differential operator of the form 

F = <T V M V„ + Q - m 2 , (1) 

acting on the (super)field <fi A (x), where 

Q = Q a b (2) 

is some matrix- valued function playing a role of the potential, V M is the appropriate covariant 
derivative and m 2 is a matrix satisfying V M m 2 and commuting with Q. It is unnecessary 
at this stage to know the exact form of the affine connection; all that is needed now is the 
knowledge of the commutator of the covariant derivatives which defines curvature according 
to a rule 

[V„,V„]<^ = lZ A B ^<f> B . (3) 

The renormalized effective action constructed from the Green function of the differential 
operator (1) is given by 

where [a n ] is the coincidence limit of the n-th HDSM coefficient and Tr is the matrix super- 
trace defined as [24] 

Tr = tr+ - tr_, (5) 

where 
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tr ± / = £/A4[l±(-irb (6) 

A Z 

and €a is the Grassman parity of c/A The coefficients a , ai, and a 2 contribute to the divergent 
part of the action , Wdiv, which have to be absorbed into the classical gravitational action 

S g = J g 1/2 d 4 x (A + X R + A XJ R 2 + X 2 R^R IX1/ ) + A 3 x, (7) 

where 

x = J g ^d A x (R^R^ - AR^RT + R 2 ) , (8) 

by renormalization of the bare coupling constant. The parameters lambda should be de- 
termined experimentally and are expected to be small, since otherwise they could cause 
observable deviations from the predictions of the general relativity. Instead of writing out 
the squared Riemann tensor we used in (7) the Gauss-Bonnet invariant, which has zero 
functional derivative with respect to the metric tensor. 

The construction of the coincidence limits of the HDSM coefficients is, except the first 
two, an extremely laborious task. The third coefficient, [a 2 ], that is proportional to the 
anomalous trace of the renormalized stress-energy tensor of the quantized massless and 
conformally invariant fields has been calculated by DeWitt [12]. The coincidence limit of 
a 3 has been obtained by Gilkey [19,20] whereas the coefficient a 4 has been calculated by 
Avramidi [21-23,25] and by Amsterdamski et al. [26]. Since we are interested in the lowest 
order of the effective action (4) we need simple and general expression for [a 3 ]. Here we use 
[03] as proposed by Avramidi [21-23] but with a different normalization: 

M 4 {p> + \{P,z m ] + \ (v„p + ij„) (w - \r) + iz w }, (9) 

where 

J, = V CT ft%, (10) 

P = Q - \\R, (11) 



5 



z { 2) = dq+ l -n pv n^ - i (^R^RT - ^r^rt^ -\°r), (12) 



and 



- R^n^n pa - i {-^° 2 R - \r^ p v u r + ^rtur^ - ^V^v^"" 

U R PR P + —R m v TtP u - -R RP ,R a P uX + J^-R ^i? A TR, ^ CT 

•j^gg P P " ~ 63 P v Q ±L pa ±L fj.uX 11 ' i 189 °" P u *1 

+ ^v^v^v.) • (13) 

In the above formulas 1 is the unit matrix, { , } is the anticommutator and we have omitted 
the field indices. 

Inserting (8) into (4) integrating by parts and making use of the elementary properties 
of the Riemann tensor one has 

W ren = 1Q2 * 2m2 J d'xg^Tr jp 3 + (R^R^ - R^ v + OR) + l -PTl^ v 
+ 1 



2 10 p 30 

1 -ROR + -Lp^uBT + -1- (-64it£P^ + 48P^P a/3 P Q /, 



630 140 p 7560 

+GR»»R% 1 R ua ^ + 17R^ aP R a f p R ap r - 28R^ p R a / p Ry i/ )} } . (14) 

This first order renormalized effective action applies to any spacetime and to any differential 
operator of the form (1). In what follows we shall confine ourselves to the operators 

+ + m 2 )0 (o) = 0, (15) 

( 7 M V M + m)0 (1/2) = 0, (16) 

(<$£□ - V„V - R£ - 5> 2 )0 (1) = 0, (17) 
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acting on the scalar, spinor, and vector fields, respectively, where £ is the coupling constant, 
and 7 M are the Dirac matrices obeying standard relations 7°7 /3 + 7^7° = 21g a/3 , and assume 
that the fields are neutral. Although neither (16) nor (17) has the form that allows direct 
application of the Schwinger-DeWitt technique the procedures described in Refs [3,24] may 
be used in this context. Specifically, by appropriate redefinition of the massive spinor field 
one obtains 



1 



m 2 ) <f>^ 



0, 



(18) 



whereas the method presented in Refs. [24] shows that the effective action of the massive 
vector field is equal to the effective action of the diagonal operator 



(<$£□ - K- 5> 2 )0« = 0, 



(19) 



minus the effective action of the massive scalar field with the minimal curvature coupling. 
The first order of the effective action is therefore 



1 



327T 2 m 2 



g^d'x < 



(0)! 



[4 

-m4 1/2) ] 

tr[aU] - [a<°> 



3|?=0J 



where the definition of the matrix supertrace has been used. 
For fields obeying Eqs. (15-17) the curvature has the form 



K 



/IV 







i-f-fR 



4 
R" 



afiu 



(20) 



(21) 



whereas inspection of (15), (17) and (18) shows that the potential matrix is 



Q 



-fiR 



(22) 



-R a , 
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Inserting (21) and (22) into (14) making use of elementary properties of the Dirac matrices 
and Riemann tensor, after simple calculations one obtains the first term of the asymptotic 
expansion of the effective action in the form [21,23] 

W ren = 7^2 2 / R^R + ^ R^UR^ + «« R 3 + Ot^ RR^ (23) 

+ ai s) RR^ pa R^ + a^R^R? + a^R^R pu R p ° v 

+ o^R^R^R"^ + a^R.rR^Rxr + a^Bf^^^ 
1 io 



192vr 2 m 2 . . 

where the numerical coefficients depending on the spin of the field are given in a Table 
I. Note that because of (20) our coefficients a^ 2 ^ for the spinor field are twice these of 
Avramidi's [^TJ, and to obtain the result for the massive neutral spinor field on has to 
multiply W^ll by the factor 1/2. 



III. THE STRESS-ENERGY TENSOR IN R = GEOMETRIES 

The renormalized stress-energy tensor is given by 
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W£ n = (T^ren, (24) 
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and for the massive scalar, spinor, and vactor fields may be rewritten in terms of the varia- 
tional derivatives of the actions Wi in the form 

Functionally differentiating the renormalized effective action with respect to the metric 
tensor, performing elementary simplifications and finally retaining in the result only the 
terms that are nonzero for R = geometries, after rather long calculations, one has 

1 " W 1 = (...), (26) 



s 
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7 — W 2 = V^R pa V l 'R (XT + V^R pa V a R pv - 3V p R pa V a R pu 

+ 2V p V u aR p p - n 2 R» u - ^V x R pa \/ x R pa g^ - V p V a nR pa g^ u 
+ 3V a V u RfR pa - V a V^R p u R pa - DR p v R p » - 2,V a V^R p(T R pu 
- URfR pv + V p V u R pa R ap + (...), (27) 

' d W 3 = (...), (28) 



9 1/2 5fiV 
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W 4 = 2V»R pa V v R pa - 2V x R P aV x R pa g pv + 2V v V p R pa R p(y 



- 2UR pu RP (T g^ - R pry R pa R pv + (...), (29) 

4jt^5 = 2V^ PCTA7 V^^ - 2V T R paXl V T R p ^g^ 
9 1/2 5g pu 

+ 2V^V^ pct a 7J R p,jA7 - 2UR paXl R p ^g^ - R^R paXl R p ^ + (...), (30) 
! ' -W 6 = ZV U R pa V a R pp - 3V a R p 1/ V 1/ R a > 1 - -V \R pa V° R px g pv 



g 1/2 5g^» p p 2 A ^ 

+ 3V a V v R p R pa - ^-V a V x R pa R px g pu + 3V CT V" R pc R m - \uR»R pp 



- ZR p °R:R pp - ^URfBr + ^R pa R x p R° x , (31) 
ttttt — W 7 = V a R p V p R au - 2V v R pa V x R x pap + 2V \R p(T V v R pXap 



g l/2 Sg»» 



- 2V x R pa V x R ppau + 2V Aj R p(T V 7 i?/ CTA - V a V p R pv R pa + 2V X V V R p p aX R pa 

- nR p ^R p ° - V a V^R paX yR px g pu + Vv^/iT" + Jvv^/F" 

+ ^i? p(TJ R A7 ^ A ^ + 2V A V^ p(T ^ A ^ - 3R pa RfR px ° u - x R pfJ R^ x g^ 

- DR pa Rp^\ (32) 



-W 8 = -2V a R p V x R x pav + V a R p V x R x pa 



9 1 ' 2 Sg^ 

+ V v R paXl V^ R pcjXp - V-yR p jV~<R paXp + V^A-yV"/^ 7 



- 2V x V p R p l, a p R uX - 2V p V x R paX v R^ + V u V 1 R/ l aX R pl ' TX 
+ RfR aX »R p ^ x - \v T V°R paXl R pTX ig^ - ±nR p v aX R"" x 
+ V°V u R paXl R p ^ - l -UR p \ x R p ™ x + 2R pa R pp R^ v 

\v r V p R paXl R^g p » - \. 
+ V x V a RfR pXc \ (33) 



+ -V T V p R paXl R° TX tg pv - -R pa R Xl p R^ - 2V p V x R pa R^ Xv 



J_ JL^ 9 = -QV x R p JV^R^ - QV,R p JV x R p ^ 
9 og pv 

- ZR p jR 1T Xv R p(J ~ <r - 2V A V 7 i^i? pCT7At - 4V x V^R p jR pa ^ 

- AV^ p R p p aX R^ - 2V,V p R p p aX R^ + X -R p(yXl R T p ° R^ K gT, 

(34) 



and 



4hY~ W ^ = WR^VrRf™ + 3V,R p JV°R p ^ 
9 ^9/j.i/ 



+ 3V p R paXl V~ / R CTfMXv + 3V p R paXl V^R auXp - 3V 7 V A i?//i? p7,TA 

+ V 1 V x R p u aX R p ~ /ap + 2V 7 V A i?/ CTA i? p7 ^ + 3R paXl R T a ' <p R pTXv 

+ 2V <T V 7 i? p<TA i/ i? wA7 - 3V CT V 7 i? pCTA7 i?^ A " + V a V 7 R p jR puXy 

- l -R paXl R T P :R^g p \ (35) 

where the ellipsis denote omitted terms that contain the scalar curvature and its covariant 
derivatives. For nondiagonal R = metrics obtained {T p )^ r s } n must be symmetrized. As 
expected, the functional derivatives W\ and W3 do not contribute to the stress-energy tensor 
in the Reissner-Nordstrom spacetime. It should be noted that Eqs ( 26-35 ) have been 
obtained by putting R = in the general result, which is more complex and shall not be 
presented here. Inspection of Eqs. (26-35) shows that to construct the stress-energy tensor 
of the massive fields in the Ricci-flat geometry it suffices to analyse only W 5 , W s , W 9 , W 10: 
and therefore our (T£)( s J n generalizes earlier results derived by Frolov and Zel'nikov [1-3]. 
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IV. {Ttf)^ ] EN IN THE REISSNER-NORDSTROM SPACETIME 



Now we are ready to construct the stress-energy tensor of the massive quantized fields 
in the nonextremal Reissner-Nordstrom geometry described by the line element 



ds " = ~ (l - ™ + ^) dt2 + (l " ™ + ^) 1 dr ' + r2 ( sin2 ed <f + d ° 2 ) > ( 36 ) 

where M is the mass and e is a charge of the black hole. For e 2 < M 2 the equation g 00 = 
has two positive roots 

r± = M ± (M 2 - e 2 ) 1 / 2 , (37) 

and the larger root represents the location of the event horizon whereas r_ is the inner 
horizon. In the limit e 2 = M 2 horizons merge at r = M, and the the Riessner-Nordstrom 
solution degenerates to the extremal one with the line element given by 

ds 2 = (l - yfdt 2 + (l - 2 dr 2 + r 2 (sin 2 ^0 2 + d6 2 ) . (38) 

Although the stress-energy tensor could be evaluated in the nonstatic background (pro- 
vided the changes of the geometry are slow) we shall confine ourselves to the exterior region 
where the spacetime is static. We begin with the massive scalar field extensively studied in 
the Ref.[8]. Constructing components of the Riemann tensor, its contractions and covariant 
derivatives and subsequently inserting them with appropriate coefficients af^ into (25) we 
have 

(T^l = CI + - I) DS, (39) 

where 

C\ = \ - - ( 1248 e 6 - 810 r 4 e 2 + 855 M 2 r 4 + 202 r 2 e 4 

* 30240 ir 2 m 2 r 12 V 

- 1878 M 3 r 3 + 1152 Mr 3 e 2 + 2307 M 2 r 2 e 2 - 3084 rMe 4 ) , (40) 

D\ = } - - f— 792 M 3 r 3 + 360 M 2 r 4 + 2604 M 2 r 2 e 2 

720 n 2 m 2 r u v 

- 1008 Mr 3 e 2 - 2712 rMe 4 + 819 e 6 + 728 r 2 e 4 ) , (41) 
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CI = -— - \ „ 10 f 444 e 6 - 1488 Mr 3 e 2 + 162 rV + 842 r 2 e 4 - 1932 rMe 



30240 7r 2 m 2 r 12 

+ 315 M 2 r 4 + 2127 M 2 r 2 e 2 - 462 M 3 r 3 ) , (42) 
D r r = } - - (-792 M 3 r 3 + 360 M 2 r 4 + 2604 M 2 rV - 1008 Mr 3 e 2 



720vr 2 m 2 r 12 

- 2712 rMe 4 + 819 e 6 + 728 r 2 e 4 ) , (43) 

°o = ~ w^n \ 2 12 (3044 rV - 2202 M 3 r 3 - 10356 rMe 4 
30240 71^ m z r lz v 

+ 3066 e 6 - 4884 r 3 Me 2 + 9909 r 2 M 2 e 2 + 945 M 2 r 4 + 486 r 4 e 2 ) , (44) 



and 



Dn = - - f 3276 r 2 M 2 e 2 -1 176 r 3 Me 2 - 3408 rMe 4 + 1053 e 6 

720 7r 2 m 2 r 12 v 

- 1008 M 3 r 3 + 432 M 2 r 4 + 832 rV) . (45) 

Obtained result for nonvanishing components of the stress-energy tensor coincides with the 
(Tjf)ren constructed by Anderson, Hiscock and Samuel. This coincidence is, of course, not 
surprising as there is a one-to-one correspondence between the order of the WKB approx- 
imation and the order of the Schwinger-DeWitt expansion. To obtain the m -2 — terms one 
has to use 6-th order WKB approximation of the mode functions and the results (39-45) are 
simply manifestation of this correspondence. 

Having computed functional derivatives of Wi with respect to the metric tensor the 
construction of the stress-energy tensor of the massive fields of higher spins present no 
problems. Indeed, inserting coefficients a^ 2 ^ for the neutral spinor field into (25) one 
obtains 

(T*)Wp = \ - 12 (2384 M 3 r 3 + 10544 r 2 e 4 - 22464 r 3 Me 2 + 21832 r 2 M 2 e 2 

- 1080 M 2 r 4 - 21496 rMe 4 + 4917 e 6 + 5400 r 4 e 2 ) , (46) 



(T;) - } = 403207T 2 m 2 r^ 2 (504 M 2 r 4 + 1080 r 4 e 2 - 784 M 3 r 3 - 6336 r 3 Me 2 

+ 3560 r 2 e 4 + 8440 r 2 M 2 e 2 - 8680 rMe 4 + 2253 e 6 ) , (47) 
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and 

(7?)£/ 2 ) = \ - - f— 3536 M 3 r 3 + 12080 r 2 e 4 - 20016 r 3 Me 2 + 30808 r 2 M 2 e 2 

40320 7r 2 m 2 r 12 v 

+ 1512 M 2 r 4 - 33984 rMe 4 + 9933 e 6 + 3240 r 4 e 2 ) . (48) 

Similarly, repeating the steps for the massive vector field one gets 

(T/)^ = \ - - (-31057 e 6 - 1665 M 2 r 4 - 41854 r 2 e 4 - 93537 r 2 e 2 M 2 

10080 ix 2 m 2 r 12 v 

+ 107516 rMe 4 + 3666 M 3 r 3 + 69024 r 3 e 2 M - 12150 e 2 r 4 ) , (49) 



(T r )W = - - f 1050 M 3 r 3 - 693 M 2 r 4 + 12907r 2 e 2 M 2 - 10448 r 3 e 2 M 

10080 n 2 m 2 r 12 v 

- 16996 rMe 4 + 2430 e 2 r 4 + 6442 r 2 e 4 + 5365 e 6 ) , (50) 

and 

{Tg)<U n = - innon \ - - ( 13979 e 6 - 2079 M 2 r 4 + 20908 r 2 e 4 + 30881 r 2 e 2 M 2 
10080 n 2 m 2 r 12 v - 

- 44068 rMe 4 + 4854 m 3 r 3 - 31708 r 3 e 2 M + 7290 e 2 r 4 ) . (51) 

Simple calculations show that the tensors (46-48) and (49-51) are covariantly conserved. 
Moreover, it could be easily verified that taking the limit e = results, as expected, in the 
formulas derived by Frolov and Zel'nikov in the Schwarzschild spacetime (see for example 
Ref. [5] ). Although there are no numeric calculations of the stress-energy tensor of the 
massive spinor and vector fields against which one could test the results (46-51), we expect 
that the approximation is reasonable so long the mass of the field is sufficiently large. 

Since the Schwinger-DeWitt approximation is local and the geometry at the event horizon 
is regular, one expects that the stress-energy tensor is also regular there. Indeed, it could 
be easily shown that if there are no fluxes of energy the regularity conditions on the event 
horizon [27,28] require that the components of the stress-energy tensor and 

2M e 2 ^" 1/2 



1 " — + 3 (TtYl ~ mfl) (52) 
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remain finite as r — > r+. Since the difference between the time and radial components of the 
stress-energy tensor factors, i.e. 

(T?>S. - = (l - ^ + j£) ^V), (53) 

where F^(r) for each spin of the field is a simple polynomial in 1/r, one can draw a general 
conclusion that our approximate the stress-energy tensors are regular as one approaches the 
event horizon. Analyses carried out in Ref . [8] indicate that all components of the numerically 
evaluated stress-energy tensor of the massive scalar field are also finite on the event horizon. 
Repeating calculations in the spacetime of the extreme Reissner-Nordstrom black hole, one 
obtains the stress-energy tensor that is regular on the event horizon. Simple calculations 
show that 

1 " -) ( 54 ) 

is finite at r = M. 

To study obtained {T")^) n further it is useful to introduce new coordinate x = (r— r + )/M 
and a new parameter q = \e\/M. Since for the massive scalar field Anderson, Hiscock, and 
Samuel have found that for m > 2/M the Schwinger-DeWitt approximation is rather good 
near the event horizon we also take this bound in our calculations of the vacuum polarization 
of the spinor and vector fields. Our results for q = 0, q = 0.95, and m M = 2 are displayed 
in the Figures 1-6. Inspection of the figures shows that for q close to the extremal value, the 
radial dependence of the components of the stress-energy tensor of the massive spinor field 
and their vector counterparts is qualitatively similar. On the other hand, for small q only 
the radial components exhibits such a similarity. Indeed, on the event horizon {T£)$ n and 
(^?)ren 2 ' ) differ in sign. The difference in the sign of the horizon values of the stress-energy 
tensor occurs also for the tangential components. Moreover, in the vicinity of the event 
horizon the magnitude of the vacuum polariation effects increases with spin of the quantized 
field. 

Geometries that could be obtained from nonextremal black holes taking the extremality 
limit and expanding the near-horizon geometry into the whole manifold received recently 
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much attention. Near the event horizon of the extremal Reissner-Nordstrom black hole the 
geometry approaches that of the Bertotti-Robinson [29] 

M 2 



^_ (_<ft2 + d ~2 + f 2 d0 2 + 



~2 • 2 

r sin 



(55) 



as could be demonstrated [28] expanding the line element (55) in power series about the 
event horizon and subsequently making substitution 



r / 



(56) 



The stress-energy tensor of the massive fields in the Bertotti-Robinson may be easily obtained 
either by constructing the curvature terms for the line element (55) and inserting them into 
(25) or taking |e| = M limit in the (T^)^ s J n near the event horizon. Simple calculations 
give |30| 



ren 



2880vr 2 m 2 M 6 



diag[l, 1,-1,-1], 



(57) 



where 



(») 



21 

37 
14 

114 
7 



(58) 



Assuming that the renormalized cosmological constant, A ren , is zero in the analog of 
the gravitational action (7) with the renormalized bare lambda coefficients, the Bertotti- 
Robinson geometry is a self-consistent solution of the semiclassical Einstein equations with 
the source term given by the stress-energy tensor of the massive field in the large mass 
limit if /i^ < 0. It is because 



9 1/2 Sg^ J 



(59) 



and 



JfJ,V 



_l 6_ 



(60) 
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vanish for the line element (55). An interesting consequence of (57) is that a self-consistent 
solution is possible for the massive scalar field provided £ > 5/14, whereas for massive 
spinor and vector fields appropriate solutions do not exist. It should be noted however, that 
the stress-energy tensor of the massive scalar field with the physically most plausible values 
of the coupling constant, namely £ = and £ = 1/6, do not yield self-consistent solutions 
and therefore the scalar field is not different than the spinor or vector field in this regard. 

V. CONCLUDING REMARKS 

In this paper we have constructed the renormalized stress-energy tensor of the massive 
scalar, spinor, and vector fields in the Reissner- Nordstrom spacetime. The method employed 
here is based on observation that the first order effective action could be expressed in terms 
of the traced coincidence limit of the coefficient a 3 . The general (T^)^ n , that has been 
obtained by functional differentiation of the effective action with respect to a metric tensor, 
consists of over one hundred terms, such as the terms cubic in curvature or involving fourth 
derivatives. Since even after simplifications, the final result is rather complicated the specific 
calculations are long but straightforward. 

Applying Eqs. (26-35) to the massive scalar field we rederived the results of Anderson, 
Hiscock, and Samuel. Their calculations were based on the WKB approximation of the 
solutions of the scalar field equation and summation thus obtained mode functions by means 
of the Abel-Plana formula. On the other hand, the method employed here may be regarded 
as geometrical and the identity of results is, although expected, impressive. To our knowledge 
spinor and vector fields have not been discussed earlier. 

The results (39-51) have also been used to construct and analyse stress-energy tensor 
in the two interesting limiting cases that could be obtained from the Reissner-Nordstrom 
solution: the extremal Reissner-Nordstrom and Bertotti-Robinson geometries. Because of 
the form of the stress-energy tensor and the fact that the variational derivatives of the 
functionals constructed from R^R^ and R 2 vanish in the Bertotti-Robinson spacetime, 
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this geometry may be a self-consistent solution of the semiclassical Einstein field equations. 
We found that the self-consistent solution is possible for the massive scalar field provided 
£ > 5/14, whereas for massive spinor and vector fields such solutions do not exist. 

Finally, we remark that it would be interesting to construct the next order of the renor- 
malized effective action (4). As the functional W ren at that order involves coincidence limit 
of the 04 coefficient, which is, in turn, given by a very complicated formula, one expect that 
such a calculation would be a real challenge. Another important direction of investigation 
is generalization of the obtained results to the elliptic operators (1) with other physically 
interesting matrix potentials Q and curvatures 1Z A b, and to analyse the back reaction of the 
quantized massive fields on the metric. We hope that the results obtained in this paper will 
be of use in further calculations. 

It should be emphasized however, that being local in its nature, the Schwinger-DeWitt 
expansion does not describe particle creation which is a nonperturbative and nonlocal phe- 
nomenon. Moreover, in the massless limit the method breaks down. To address successfully 
this group of problems new methods are necessary as, for example, the covariant perturba- 
tion theory [ ]5TJ . 
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TABLES 



(s) 

TABLE I. The coefficients aq for the massive scalar, spinor, and vector field 
s = s = 1/2 s = 1 

lc2 _ U , 1 3_ 27 

2^ 5^ ' 56 140 280 

1 J_ _9_ 

140 14 28 



1 __5_ 

432 72 



30 V 6 W 90 60 



30 V 6 

J-fl-^ __ Z_ — -L 

30 V 6 > J 720 10 

8 25 _52 

945 376 63 

2 47 19 

315 630 105 

1 19 61 

1260 630 140 

17 29 67 

7560 3780 2520 

1 _ J_ J_ 

270 54 18 



18 



REFERENCES 

[1] V. P. Frolov and A. I. Zel'nikov, Phys. Lett. B 115, 372 (1982) 
[2] V. P. Frolov and A. I. Zel'nikov, Phys. Lett. B 123, 197 (1983) 
[3] V. P. Frolov and A. I. Zel'nikov, Phys. Rev. D 29, 1057 (1984) 

[4] V. P. Frolov, Proceedings of the Lebedev Institute of the Academy of Science of USSR 
169 (1986) 

[5] V. P. Frolov and I. D. Novikov, Black Hole Physics (Kluwer Dordrecht 1998) 
[6] L. A. Kofman and V. Sahni, Phys. Lett. B 127, 197 (1983) 
[7] P. R. Anderson, Phys. Rev. D 41, 1152 (1990) 

[8] P. R. Anderson, W. A. Hiscock, and D. A. Samuel, Phys. Rev. D 51, 4337 (1995) 
[9] W. A. Hiscock, S. L. Larson, and P. R. Anderson, Phys. Rev. D 56, 3571 (1997) 

[10] B. E. Taylor, W. A. Hiscock, and P. R. Anderson, Phys. Rev. D 55, 6116 (1997) 

[11] R. Herman, Phys. Rev. D 58, 084028 (1998) 

[12] B. S. DeWitt, Dynamical theory of groups and fields (Gordon and Breach, New York 
1965) 

[13] S. M. Christensen, Phys. Rev. D 14, 2490 (1976) 
[14] S. M. Christensen, Phys. Rev. D 17, 946 (1978) 
[15] B. S. DeWitt, Phys. Reports 19 C, 297 (1975) 
[16] J. S. Schwinger, Phys. Rev. 82, 664 (1951) 

[17] Birrell and P. C. W. Davies, Quantum fields in curved space (Cambridge University 
Press, Cambridge 1984) 



19 



[18] S. A. Fulling, Aspect of quantum field theory in curved space-time (Cambridge University 
Press, Cambridge 1989) 



[19 
[20 
[21 
[22 
[23 
[24 
[25 
[26 

[27 
[28 
[29 

[30; 

[31 



P. B. Gilkey, J. Diff. Geom, 10, 601 (1975) 

P. B. Gilkey, Trans. Am. Math. Soc. 225, 341 (1977) 

I. G. Avramidi, |hep-th/9510140| 



I. G. Avramidi, Teor. Mat. Fiz, 79, 219 (1989) 
I. G. Avramidi, Nucl. Phys. B 355, 712 (1991) 

A. O. Barvinsky and G. A. Vilkovisky, Phys. Reports 119, 1 (1985) 
I. G. Avramidi, Phys. Lett. B 238, 92 (1990) 

P. Amsterdamski, A. L. Berkin, and D. J. O'Connor, Class. Quantum Grav. 6, 1981 
(1989) 

C.-G. Huang, Phys. Lett. A 164, 384 (1992) 

P. R. Anderson, W. A. Hiscock, and D. J. Loranz, Phys. Rev. Lett. 74, 4365 (1995) 

B. Bertotti, Phys. Rev. 116, 1331 (1959); I. Robinson, Bull. Acad. Pol. Sci. 7, 351 
(1959) 

Our /i < -°- ) differs from the value evaluated in Ref. which is insenstive to the coupling 
constant £. 

G. A. Vilkovisky, in Quantum Theory of Gravity, edited by S. M. Christensen (Hilger, 
Bristol, 1984). 



20 



FIGURES 




FIG. 1. This graph shows the radial dependence of the rescaled component (Tf)ren^ 
(A = 180(8M) 4 7r 2 ) of the renormalized stress-energy tensor of the massive spinor field with 
m = 2/M. From top to bottom the cureves are for q = 0.95 and q = 0. 
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FIG. 2. This graph shows the radial dependence of the rescaled component (TJT^en^ 
(A = 180(8M) 4 7r 2 ) of the renormalized stress-energy tensor of the massive spinor field with 
m = 2/M. From top to bottom the cureves are for q = 0.95 and q = 0. 
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FIG. 3. This graph shows the radial dependence of the rescaled component (Tfj)^^ 
(A = 180(8M) 4 7r 2 ) of the renormalized stress-energy tensor of the massive spinor field with 
m = 2/M. From top to bottom the cureves are for q = and q = 0.95 
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FIG. 4. This graph shows the radial dependence of the rescaled component (Tjf)ren 
(A = 90(8M) 4 7r 2 ) of the renormalized stress-energy tensor of the massive vector field with 
m = 2/M. From top to bottom the cureves are for q = 0.95 and q = 0. 
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FIG. 5. This graph shows the radial dependence of the rescaled component (T£)£) n 
(A = 90(8M) 4 7r 2 ) of the renormalized stress-energy tensor of the massive vector field with 
m = 2/M. From top to bottom the cureves are for q = 0.95 and q = 0. 
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FIG. 6. This graph shows the radial dependence of the rescaled component (Tg)rJ n 
(A = 90(8M) 4 7r 2 ) of the renormalized stress-energy tensor of the massive vector field with 
m = 2/M. From top to bottom the cureves are for q = and q = 0.95 
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